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As a test of the standard inflationary cosmology, which generically predicts nearly scale-invariant 
spectrum of primordial curvature fluctuations, we perform Markov-Chain Monte-Carlo analysis to 
search for possible modulations in the power spectrum and determine its shape together with the 
cosmological parameters using cosmic microwave background radiation data. By incorporating var- 
^ ' ious three-parameter features on the simple power-law spectrum, we find an oscillatory modulation 

localized around the comoving wavenumber k ~ 0.009Mpc _1 at 99.995% confidence level which im- 
. proves the log-likelihood as much as — A21n£ = AxcB = —22. This feature can be detected even if 

■ we use only the cross correlation between the temperature and the E-mode polarization anisotropies. 

(N 



PACS numbers: 98.70.Vc, 95.30.-k, 98.80.Es 



I. INTRODUCTION 
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' In the standard inflationary cosmology p], E|, H, Q the observed large-scale structures and the anisotropics in 

the cosmic microwave background radiation (CMB) originate in tiny quantum fluctuations generated during the 
inflationary expansion stage in the very early Universe 0, @, 0| ■ Conventionally, the power spectrum of primordial 
curvature fluctuations, P(k) = (\lZk\ 2 ), has been assumed to follow a simple functional form such as a power-law 
A(k) = -P(fc) = A{k/ko) n3 ~ 1 which is quantified by the amplitude A and the spectral index n s , and the previous 
statistical analyses of the power spectrum have mostly focused on these two parameters (and at best the scale 
dependence of n s , the running, which may be important to distinguish the inflation model[8|). It is true that the 
simplest class of inflation models predicts a power-law spectrum with n s close to unity 0, 0,0] , but we may not be 
able to identify the correct theory of the early Universe if we restrict our parameter space from the beginning. From 
the viewpoint of observational cosmology, the shape of the primordial fluctuation spectrum should be determined 
purely from the observational data without any theoretical prejudices. 

Along this line of thought, several Markov-Chain Monte-Carlo (MCMC) analyses have been performed to search 
for possible deviation from the power law, but none has detected statistically much significant modulations so far 
EOj EH, Ell [H, 0] • To be more quantitative, the presence of a nontrivial feature may be decided in terms of Akaike's 
Information Criterion (AIC) (l5l | which asserts that introduction of an additional fitting parameter is justified if and 
only if Xefi improves more than —2 with it. Most of the previous analyses resulted in the improvement of ^xts mucn 
less than that required by AIC, and others reached it only at a marginal level, so that it was hardly possible to claim 
■ the presence of a feature. It should be emphasized, however, that those analyses did not have sufficient resolution to 
detect spectral fine structure, because they tried to fit the primordial power spectrum in a broad range of scales in 
terms of a limited number of degrees of freedom using sparse sampling or some specific theoretical models. This lack 
of resolution was inevitable because exploring large dimensional parameter space is extremely time-consuming even 
with the help of MCMC analysis. 

Recently, on the other hand, two of the present authors (RN and JY) performed reconstruction of primordial power 
spectrum from the angul ar p ower spectrum of the CMB temperature fluctuations, Cj T , of the five-year WMAP 
data (hereafter WMAP5) [l6l. [l7j using two different non-parametric methods, namely, the cosmic inversion method 
0, EH, E3, and the maximum- likelihood reconstruction method 0, 0, E||. They have probed the 
primordial power spectrum with finer resolution than aforementioned MCMC analyses and found an anomalously 
large deviation from the best-fit power-law primordial spectrum of WMAP5 around the wavenumber k ~ 0.009Mpc _1 
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or kd ~ 124 where d = 1.43 x 10 4 Mpc is the distance to the last scattering surface. This scale corresponds to the 
multipole £ ~ 120 and the length scale r ~ 710Mpc. 

In their reconstruction procedure, however, they had to fix values of the cosmological parameters and adopt an ad 
hoc smoothing prior to ensure sensible reconstruction. The latter tends to discard a large part of the information on 
the fine structure which may have affected the significance of the claimed spectral feature itself. Furthermore although 
the reconstructed power spectrum exhibited oscillatory features with both peaks and dips, the result of decomposition 
of the continuous reconstructed curve to statistically independent band-powers, which is necessary to discuss their 
statistical significance, revealed only a 3.3<r peak but no dips (26[. 

The purpose of this paper is to focus on the fine structures of the spectral shape. We use a forward analysis 
implemented by MCMC simulation which circumvents the above-mentioned difficulties peculiar to the inverse mapping 
approach. We also overcome the lack of resolution, which was a very serious problem of previous MCMC analyses 
as mentioned above, by searching for and concentrating on the most prominent localized feature imprinted on the 
otherwise power-law spectrum. 

II. METHOD 

First we prepare various model power spectra in which a localized spectral feature is imprinted on the power-law 
with three fitting parameters. We then calculate the temperature-temperature (TT) angular power spectrum and 
the cross correlation between the temperature and the E-mode polarization (TE) to find the best-fit values of these 
parameters as well as those of other cosmological parameters in the ACDM model using the CosmoMC[27| code. 

We start with the following three types of three-parameter modulations. 

A-type: a A-shaped peak is introduced at k = fc* and connected to a power-law A(k) = A(fc/fco)" s_1 by straight lines 
in the logfc — A(k) plane. The location of the peak, k*, its width, and height are additional fitting parameters. 
Here ko = 0.002Mpc _1 is the pivot scale. 

\A-type: a peak at k — k* and a dip at k < fc* with the same amplitude of deviation from the standard power-law 
are connected by picccwise straight lines to A(k) = A(k/ko) ns ~ 1 in the logfc — A{k) plane. Again the height 
and the width of the peak are the parameters characterizing modulation besides the location k* . 

5-type: peaks and dips are characterized by the following smooth function. 

A(k) = Aik/ko) 71 *- 1 + B{k/k ) n '- 1 exp[-(fc - K) 2 /k 2 } cos[7r(fc - fe»)//c], (1) 

where A?*, B, and k are additional fitting parameters. 

We assume flat priors for the additional parameters to cover the parameter range as broad as possible. For the 5-type 
modulation, for example, we take 10 10 B = [0, 10 3 ] and 10 4 k = [l,10]Mpc _1 . As for the locaiton of the feature, fc», 
we started our calculation allowing it to vary in the full range of obscrvationally accessible domain with a sufficient 
accuracy so that we could detect a modulation localized anywhere in k space. We have confirmed, however, that in 
the range of the wavenumber accessible by current observation the most prominent feature is located at k*d ~ 124 as 
was found in the reconstruction analysis [26[ together with other possible modulations whose statistical significance 
is smaller. In fact, if we allowed fe* to vary beyond k*d > 135 and/or k*d < 115, other features would also contribute 
to MCMC analysis, albeit with little significance, in addition to the most prominent one at k*d ~ 124. Then our 
three-parameter models would no longer be a good parametrization and the MCMC calculation would not converge 
properly. We therefore decided to limit the range of fc* to k t d = [115, 135] in the final simulations in order to calculate 
the probability distribution of the amplitude of the most prominent modulation, B, properly. 

In order to check the stability of our result with respect to other features localized at different wavenumbers we ran 
MCMC analysis incorporating another feature beyond k*d > 135 at the same time with three more fitting parameters. 
As a result we found that the inclusion of additional parameters to incorporate another feature does not affect the 
posterior distributions of the original three parameters for the most prominent feature at k*d ~ 124, which implies 
that we can treat the relevant feature independently from other possible features at different wavenumbers. 

Based on the above observations we have also adopted another model to fit the feature around the most prominent 
modulation at k*d ~ 124 as follows. 

VF-type: three-step modulation on a power law with the functional form, 

A(k) = Aik/ko)"*- 1 + B < 6{kd - 114)61(122 - kd) 

+B 9{kd - 122)0(126 - kd) + B>6(kd - 126)0(134 - fed), (2) 

where B, B <y and B > are fitting parameters. 
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FIG. 1: Posterior distribution functions for parameters of the modulation amplitude B. Here the horizontal axis denotes the 
parameter B. For all models the simple power law spectrum corresponds to B = 0. 

Among the four models we use, the A-type modulation contains only an extra power so that it enhances the disper- 
sion of fluctuations compared with the background power-law spectrum for the same value of A. In the other three 
types of models, both excess and deficit are incorporated so that they do not necessarily affect overall normalization 
of fluctuations. In these three models, we put A(k) = whenever A(k) takes a negative value. In the course of 
MCMC calculations, we have found a degeneracy between the width and the height of modulation in A- and \A-type 
models and they tend to give a larger modulation amplitude with a smaller width. Such a tendency is undesirable 
from physical point of view because A(k) should be positive definite. Hence we fixed the width of modulation to 
A In A: — 0.043 in these two models which is the best-fit value in our first trial run and this choice is consistent with 
the result of the S'-type model where the width k is treated as a free fitting parameter. 

III. RESULTS 

The result of MCMC calculations is summarized in Table 1. First by incorporating an excess with the A-type model, 
the effective xts = — 2 ln£, where C is the likelihood function, improves by Axis = —6.5 with three extra parameters 
added to the standard power-law ACDM model. To this end this model satisfies the AIC. Inspecting the details of 
Axoff , however, we find that the improvement of —6.5 is entirely due to the better fit to the TE data. Although the 
fit to TT data improves around I w 120, it gets worse for t < 100. This is because the transfer function from P{k) to 
Cj T is non- vanishing for t < kd so that an excessive power around kd = k*d increases all Cj T for £ < 120 to affect 
overall normalization in an unwanted manner. 

Therefore even if it was not detected in a band-power analysis (2(| , we should take neighboring dips observed in the 
reconstructed curve as in the other three models, where xts improves as much as ^XeS = ~ ~^2 with the same 
numbers of extra parameters. We can regard it extremely significant and interpret that the result strongly suggests 
a sharp and strong oscillatory deviation from a power-law around k = k*. The marginalized ID distributions of the 
parameter B made by MCMC calculations are shown in figure 1. Based on the posterior statistical distribution of 
the amplitude of spectral modulation, we find that the pure power-law model with B = is 4. Oct (4.1ct) away from 
the mean value of B for 5-type (\A-type) model, respectively. 

We also examined the probability to find such a large deviation by accumulating more samples to explore the tail 
of the posterior probability distribution for the 5-type model which improved xts the best among the models we 
adopted. We have found 304 samples in the smallest bin of the modulation parameter with B < 10~ 10 ~ A/23, from 
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FIG. 2: Angular power spectrum of CMB temperature anisotropy in the range 100 < £ < 140. Data points are from WMAP5 
with smaller error bars indicating only the measurement errors. Green line is the result of best-fit power-law ACDM model, while 
the most wavy black curve and the less wavy magenta curve represent the results of S-type modulation and A-type modulation, 
respectively. The dotted curve is the angular power spectrum recalculated from P(k) reconstructed by the maximum-likelihood 
reconstruction method 12611. 



6279082 MCMC samples generated with temperature parameter T = 1 [28|. Hence the relative probability is 

P(B < 1 x 1(T 10 ) = 4.8 x 10~ 5 . (3) 

In order to check whether we have explored the tails of the distribution with sufficient accuracy, we have ran simulations 
with different temperature parameters (T = 1.5 ,2) in the MCMC analysis and confirmed that the probability converges 
to P = (3.9 x 10~ 5 ,5.6 x 10 -5 ), respectively. Therefore we conclude that the posterior probability for B = to be 
the case is only O(10 -5 ). 

We note that among the improvement of Ax^ ff = —22 for the ifj-type spectrum, —14 is due to the improvement of 
the fit to the TT power spectrum, Cj T , while the remaining —8 is from the TE cross correlation, Cj E . Figures 2 
and 3 depict Cj T and Cj E in the relevant range. It is intriguing that significant portion of the improvement of the 
likelihood comes from the TE data which was not used in the reconstruction approach at all (26| . Therefore this may 
be regarded as an independent support to our discovery of a non-power-law feature in the primordial power spectrum. 

In order to pursue it, we have performed two distinct sets of MCMC calculations one using only TT data and the 
other TE data. In both cases the 5*-type modulation was used and the cosmological parameters have been fixed to the 
best- fit values for this model shown in Table 1. To ensure the stability of calculation we discarded any A(k) reaching 
negative values in the calculation instead of setting A(k) = 0. As a result we have obtained the following values of 
the model parameters from TT and TE data, respectively. 
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where the errors represent 95% upper and lower bounds. Here the location of the peak is determined in good shape in 
both cases and the results are in good agreement with each other. On the other hand, TE data fails to constrain the 
width parameter k and its upper bound is essentially determined by the prior cutoff. Finally as for the amplitude of 
modulation B, both datasets give similar posterior distributions and we find the relative frequency to find vanishingly 
small modulation is 0(1O -3 ) in both cases. More precisely, we find 



P(B 



TT < 1 x lO" 10 ) 



= 2 x 10~ 3 and P(B TE < 1 x 10~ 10 ) = 1 x 10" 



Thus we can conclude both TT data and TE data suggest the existence of a feature at the same wavenumber separately 
and mutually consistently. 

Finally we note that such a feature also shifts the best-fit values of other cosmological parameters slightly as seen 
in Table 1. Although the shifts arc smaller than the observational errors of five-year WMAP data, some of them are 
larger than the expected errors in the Planck mission.lt has been already shown that globally non-power-law spectra 
generally induce shifts in the estimated values of the cosmological parameters compared with the standard power-law 
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FIG. 3: Angular cross correlation of CMB temperature and E-mode polarization in the same range. Data points are the 
observational result of WMAP5. Nearly straight green line is the result of best-fit power-law ACDM model, while the black 
curve and the magenta curve represent the results of S'-type modulation and A-type modulation, respectively. 
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TABLE I: Cosmological parameters estimated from CosmoMC code using four types of modulated power spectra. Here Of,, Q m 
and Q,a are fractions of cosmic energy density in baryons, matters and cosmological constant, respectively, Ho is the Hubble 
parameter in unit of km/s/Mpc, and r is the optical depth to the last scattering surface. For A- and v A -type models the 
parameter B stands for the difference between the peak amplitude and the power- law background value at k — k*. Axis 
represents improvement of goodness-of-fit xls for the best-fit model of each modulated spectrum from the most-updated best- 
fit power-law ACDM model. A max is the maximum of the difference of each parameter between the power-law model and 
models with modulated spectra except for the A-type which does not improve the fit. Column WMAP5 (Planck) represents 
errors observed (expected) by five-year WMAP (Planck). Errors expected in Planck mission are taken from the "Blue book". 
The mean value of k in the S-type model is k = 3.58 x 10~ 4 . 



ACDM model 0, [T3, E3, H^, H3] because there is a degeneracy between the global shape of P(k) and the cosmological 
parameters whose values are imprinted on the shape of the acoustic oscillation. Nevertheless it is intriguing that 
features localized in a narrow range of wavenumber as discussed here also changes their values at a non-negligible 
level for the next-generation analysis. 



IV. CONCLUSION 



In this paper we have found a spectral feature with a 4cr deviation from the best-fit power-law primordial power 
spectrum in a narrow range around k ~ 0.009Mpc _1 , whose existence is suggested separately by both TT data and 
TE data from five-year WMAP data. This feature can be further tested and hopefully confirmed by forthcoming EE 
data from the Planck mission. 

As emphasized above, the detected feature consists not only of a peak but also a dip. Since we cannot lower the 
amplitude of the power spectrum by any sources uncorrelated with the preexistent fluctuations, this means that the 
feature was created together with the bulk of the power-law spectrum, suggesting some nontrivial phenomenon during 
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inflation. At the moment, while a number of inflation models have been proposed to produce spectral modulation, we 
are not aware of any model that can realize a spectrum with such a localized feature. Hence efforts should be made 
toward construction of models that can account for the detected fine structures. 

Currently inflation models are observationally constrained in terms of a small number of discrete parameters such 
as n s , dn s /d]nk, the tensor-to-scalar ratio r, and a nonlinearity parameter Jnl- It is difficult to determine the 
underlying particle-physics model with these parameters alone, because there are huge degeneracies. If our finding 
is further confirmed, it can certainly serve as an important clue to single out the correct model. On the other hand, 
whether this feature is a result of some nontrivial physics or a realization of an extremely rare event, we must take it 
into account in cosmological parameter estimation in the analysis of forthcoming higher precision data. 
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